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We discuss insights that may be drawn from our recent 2 flavour O (a) -improved Wilson quark simulations. We 
discuss the evidence of the onset of chiral logarithms in the pion decay constant. An overview is given of current 
extrapolation methods and a modification of chiral perturbation theory is presented as an approach for sensibly 
extrapolating to the physical quark masses. 
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1. Introduction 

In a recent work we presented results from 
the lightest UKQCD Wilson quark simulations 
for the mass of the singlet pseudoscalar meson 
and the pseudoscalar (pion) decay constant. In 
this work we shall expand upon the pion decay 
constant extrapolation details and suggest a sim- 
plistic first attempt at an improved extrapolation 
method. 

The motivation for this recent calculation of 
the pseudoscalar decay constant, /pg, is the ex- 
pectation of observing behaviour predicted by 
Chiral Perturbation Theory (x?T) as the pseu- 
doscalar mass is lowered. In particular evidence 
of the onset of chiral logs in observables. The low- 
est order chiral Lagrangian has non-analyticity 
resulting from loop corrections introduced in the 
parameters. This non-analytic behaviour pro- 
vides a good check that the lattice calculation is 
in the regime where chiral perturbation theory is 
valid. 

Computational constraints have forced us to 
work at finite lattice spacing, however the choice 
was made to compare our results to the contin- 
uum predictions of xPT. Whilst Chiral perturba- 
tion theory has been formulated for finite lattice 
spacing |2I3| , it is at the cost of additional param- 
eters. 



2. Calculating /pg 

Our results for the pseudoscalar decay constant 
were extracted from UKQCD simulations at three 
different quark masses: k = 0.1358, 0.1355 and 
0.1350. All lattices had 16 3 x 32 volumes with 
(3 = 5.2 and a non-perturbatively improved clover 
action. Using quark propagators with sources on 
the time planes t = 0, 7, 15 and 23 allowed an 
improvement in the statistics of the simulation. 

L « 4 in each of the cases we 
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As the value of i 
would expect some finite volume effects in our cal- 
culations. These are taken into consideration 0] 
in the subsequent analysis. 



Table 1 

The raw lattice value of a/pg is given by using the 
order a improved expression (1 + bAm PS )(afA + 
CAO-fp) and we tabulate these two contribu- 
tions, m 



K 


afA 


afp 


0.1358 


0.0829(26) 


0.1457(78) 


0.1355 


0.1055(14) 


0.1835(44) 


0.1350 


0.1336(11) 


0.2468(33) 



* Current Address: Physics Division, Argonne National 
Laboratory, Argonne IL 60439, USA. 



Our results for /pg were presented in |T] and are 
reproduced here in Tableland Fig.^ The lattice 
value of a/ps used in Fig. ^ is found by using the 
order a improved expression (1 4- 6A^pg)(a/A + 
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cacl/p), where /a,p & r e the axial- vector and pseu- 
doscalar decay constants and their contributions 
are listed in Tabled 
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Figure 1. The pseudoscalar decay constant in 
units of ro from UKQCD and JLQCD versus k.[Q 



To aid comparison with the only available com- 
parable data, that of JLQCD jS] , we have used the 
same perturbative formulation of the corrections 
and also the same prescription for Za, the renor- 
malisation factor. Finally as there is a discrep- 
ancy in the prescription for evaluating ro we apply 
our determination of ro to their data. The com- 
parison of results is shown in Fig. ^ The agree- 
ment at the lightest common n values is pleasing. 
The exciting feature of our new results is the in- 
dication of curvature versus m q . It is this result 
that has motivated the current investigation. 

3. Extrapolation Comparisons 

The limitation of all current lattice calculation 
of the pseudoscalar decay constant to large quark 
masses necessitates some form of extrapolation to 
the physical quark (or equivalently pion) mass. 
Over the last few years it has become accepted 
within lattice groups that some form of chirally 
motivated extrapolation is appropriate |6I7I8I9| . 
by using insights from chiral perturbation theory, 
when extrapolating physical observables. The sit- 



uation is no different in the case of /ps, how- 
ever opinions of what is appropriate differ between 
groups. 

The chiral perturbation prediction for the pion 
mass dependence of /ps is |4I10| : 
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where f w (0) is the pseudoscalar (pion) decay con- 
stant in the chiral limit. 
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Figure 2. The pseudoscalar decay constant in 
units of r from UKQCD versus the squared pseu- 
doscalar meson mass. Also shown is an expression 
including chiral perturbation theory terms to or- 
der mp s which has been fitted (see ref. 4 where 
we use /i = 0.75 GeV and rp(fi) — —2) to agree 
with the experimental values of f v and fx which 
are shown (*). An estimate [2] of the finite size ef- 
fect expected from chiral perturbation theory (to 
order m,p g ) is shown by the vertical lines. 



Recent work 0] has provided a way of estimat- 
ing the expected finite size effects. The quark 
loops generate logarithmic corrections, but more- 
over they are the source of the finite size effects. 
Using L = 1.5 fm and m ff = 400 MeV, which are 
close to our values a suggestion of the size of the 



3 



corrections may be determined. The solid curve 
shown is an expression including chiral perturba- 
tion theory terms to order rop s which has been 
fitted (see Ref. 0] where we use /i = 0.75 GeV 
and t\f(/i) = —2) to agree with the experimental 
values of /„• and fx which are shown (*). An esti- 
mate J2J of the finite size effect expected from chi- 
ral perturbation theory (to order mp S ) is shown 
by the vertical lines. As an experiment, we also 
continued the xPT prediction for /ps to almost 
twice the kaon mass (dashed curve). It is clear the 
lattice simulation behaviour in this regime is fun- 
damentally different from the naive predictions of 
%PT. This discrepancy indicates xPT (as formu- 
lated) has a limited range of overlap with current 
simulations. 

JLQCD 

Recent progress in the extension of x?T has 
been attempted by the JLQCD collaboration 0- 
Their approach used xPT as motivation for their 
functional form which, in the dynamical case 
(K S e a = KVai) ma y be written as: 

/ PS = A + Bmls + Cm 4 PS (2) 

The analytic (even powers of m ps ) terms of 
Eq. are consistent with both xPT and the 
JLQCD data set, as seen in Fig. |3J The success 
of such a fit may be seen to vindicate "leaving the 
problem of the chiral logarithm [. . . ] for future 
publication." 
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Figure 3. Chiral extrapolation of the pseu- 
doscalar meson decay constant. The fit is made 
using Eq. 



The problem with this approach is the selective 
nature of what parts of xPT are used. The chiral 
logarithm term that occurs with a power of the 
pseudoscalar mass less than to ps is ignored. Yet 
this is the greatest contribution to the curvature 
in the small mass regime. The form is also in- 
consistent with the large mass behaviour of /ps 
which experience has shown to be proportional 
to mpg up to pseudoscalar masses of the order of 
1.4 GeV 2 . Thus, any perceived successes of this 
form of extrapolation must be tempered by the 
clear inconsistencies with the behaviour of QCD 
as established by xPT. 

Dyson Schwinger Equation 

An interesting alternative insight is that pro- 
moted by the Dyson-Schwinger Equation (DSE) 
community. The Dyson-Schwinger equations, 
whilst maintaining a Poincare covariant frame- 
work, provide a setting for a non-perturbative 
chiral symmetry preserving truncation that leads 
to an efficient one-parameter rainbow-ladder 
model In this DSE approach, the pion 

is well-understood and is essentially model- 
independent; it is both the Goldstone boson and 
a bound state of strongly dressed quarks |12I13| . 
This reduction of QCD to an understandable case 
benefits us as the circumstances of the lattice 
(that is equal mass quarks as discussed later) may 
be investigated. 

The DSE approach of Tandy Jl] to the pseu- 
doscalar decay constant is presented in Fig. 0] by 
the solid curve. For the dashed curve, the same 
calculated decay constant is plotted against the 
pseudoscalar mass obtained from m q via the Gell- 
Mann-Oakes-Renner relation, as would be used 
by xPT- This alternative approach differs from 
%PT and the lattice calculations as quark loops 
are explicitly excluded from the calculation. It 
would thus be expected that the exclusion will, in 
the light quark mass region, result in a behaviour 
similar to that expected by the JLQCD collab- 
oration in Eq. (|2J) — no logarithms, and yet be 
suggestive of the lattice calculations in the region 
where quark loops become suppressed. Tandy pa- 
rameterised the DSE results with the following 
phcnomcnologically motivated form 
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Figure 4. The pseudoscalar decay constant ver- 
sus the square of the pseudoscalar mass both cal- 
culated within the Dyson-Schwinger framework. 
The solid line is the exact pseudoscalar mass, 
whilst the dashed curve is the mass obtained 
by using the Gell-Mann-Oakes-Renner relation- 
ship. [TH] 
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Whilst reproducing the correct QCD behaviour at 
large mp S it is incorrect at small quark masses. 
As was discussed above this should not be consid- 
ered a failure, as the model that motivated this 
form lacks the chiral behaviour that induces the 
logarithm in Eq. QJ. 

3.1. Modified xPT 

The final approach we will mention is a mini- 
malistic modification of xPT as motivated by our 
previous work |fil7l8llfi| . In these investigations 
we found that a change in the form of the reg- 
ularisation of xPT admits an extended radius of 
convergence whilst maintaining an exact agree- 
ment with the results of the dimensional regu- 
larised solution. Recent work of the Adelaide 
group has show that this exact agreement is 
indeed obtainable between different formulations 
of the regulariser, and it is indeed a requirement 



of a regularisation scheme. 

The changes we make to the dimensional regu- 
larised result of xPT to form our Modified xPT 
are two: 

• Terms of 0(rrip S ) and higher in the full ex- 
pression of Eq. are neglected, 

• An additional parameter Bmp S is intro- 
duced to absorb the physics ignored by the 
truncation. 

The form used for extrapolating thus becomes: 
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where the new term BmL suppresses the chiral 
logarithms above a certain energy scale and pro- 
vides the missing analytic behaviour near the chi- 
ral limit. A similar approach has been attempted 
by JLQCD Ql|. 

The advantages of this form for extrapolating 
results from lattice QCD calculations are that it 
reproduces xPT near the chiral limit, but it also 
has the intermediate mass behaviour of QCD that 
the lattice has shown. This simple form allows an 
extension of the radius of convergence of xPT- It 
must be noted that this form is indeed naive in the 
way the two mass limits are enforced but it does 
suggests hope for a more rigorous FRR solution 
that has been achieved elsewhere ^7j. Such an 
extension is part of an ongoing investigation. 

This modified xPT approach is more suited 
to extrapolating lattice QCD calculations, at the 
current juncture, than the dimcnsionally regu- 
larised approach of Eq. for another reason too: 
it has fewer fit parameters. As shown in Fig. 2] 
near the chiral limit we currently have, at best, 
5 data points, and fitting purely to our UKQCD 
data we would have only three data points. The 
0(p 6 ) low-energy constants expansion of Eq. (JTJ 
(see Ref. 0] for the full expression) has seven pa- 
rameters: /tt(0), Ai, . . . , A 4 , r F (n) and n; each 
of which must be fit to some data. This has been 

FT 



undertaken in 



but in the comparison of 



lattice QCD to experimental results, these pa- 
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Figure 5. The pseudoscalar decay constant in units of ro from UKQCD versus the squared pseudoscalar 
meson mass. The experimental values of / w and fx are shown (*). The dotted curve is the fit of 
the Modified %PT result of Eq. (0J. The vertical lines indicate expected finite size effects at various 
pseudoscalar meson masses, as in Fig. |3] 



rameters must independently be determined by 
the theory. 

Conversely the simplistic approach of Eq. 
has only three parameters /tt(0), A and B. In 
comparing these two approaches the expectation 
that the non-analytic behaviour should be the 
same is an important success. The values of A 
(Eq. (0J) and A 4 (Eq. QJ) are consistent within 
the large uncertainties as previously discussed 
in 0. The value of A 4 from 4 is 1.25 to.ls GeV 
whilst we find A = 0.93 ± 0.42 GeV. In a simi- 
lar vein to Ref. some of the physics we have 
neglected, i.e. higher powers of m PS , has been 
absorbed into the parameter B. We present our 
fit of Eq. gj) to our UKQCD data in Fig. [SJ The 
vertical dashed lines are the same finite size error 
estimation calculations as discussed in Sec. |3J 

It is clear that even with a finite size correc- 
tion to the UKQCD calculations there is a dis- 
crepancy between the data obtained from lattice 
simulations and the experimentally measured de- 



cay constant of the K meson. This discrepancy 
highlights some of the limitations of the current 
results. The finite size corrections that have been 
applied to this data are indicative only and do not 
necessarily account for all the finite size effects in 
the simulation. The benefit is that they reinforce 
the intuition that finite size effects become less 
important as the quark masses get larger. 

The discussion herein of an extrapolation of the 
quark mass ignores the other important extrapo- 
lation that has not been undertaken in these sim- 
ulations; that is the continuum extrapolation. In 
this work we have chosen a value of ro as sug- 
gested by |21I22| . however the uncertainty of this 
value is 5%. Additionally, the value of Za we use 
to renormalise our result is determined perturba- 
tively. At first order the correction is 25% and 
the systematic error at second order is expected 
to be up to 5%. Whilst the value of Z^ remains 
determined perturbatively we would expect that 
part of the discrepancy in our results to remain. 
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Finally there is a fundamental difference be- 
tween the object that is calculated at non-zero 
quark masses and the pseudoscalar meson that 
is experimentally detected. The K meson con- 
ceptually is a strange quark and a light quark in 
a sea of light quarks. The pseudoscalar meson 
calculated on the lattice with the same mass con- 
sists of two quarks with masses half of that of 
the strange quark in a sea of equivalently heavy 
quarks. The DSE calculations shown in Fig. 0] 
have this behaviour, but the experimental points 
of Figs. [3 and are different objects. There is an 
expectation that the difference between a (qQ) 
object and (^^} meson is not large, but this is 
yet to be quantified. 

4. Conclusion 

We have presented a review of our recent pub- 
lication £Q and in particular the discussion of the 
extraction of the pseudoscalar decay constant on 
the lattice. We have shown that this data set 
contains an indication of the onset of chiral log- 
arithms for the first time in this quantity. Ad- 
ditionally we have discussed how the need to ex- 
trapolate to the chiral limit still exists and fur- 
thermore that the suggestion of curvature neces- 
sitates the use of chiral perturbation theory. Fi- 
nally we discuss some approaches and present a 
first attempt at a form that respects chiral per- 
turbation theory whilst being applicable in the 
region where the lattice calculations are still be- 
ing performed. 
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